Spin polarized transport driven by square voltage 
pulses in a quantum dot system 



F. M. Souza and J. A. Gomez 

International Centre for Condensed Matter Physics, Universidade de Brasilia, 
70904-910, Brasilia-DF, Brazil 

E-mail: fmsouza@unb.br 

Abstract. We calculate current, spin current and tunnel magnetoresistance (TMR) 
for a quantum dot coupled to ferromagnetic leads in the presence of a square wave 
of bias voltage. Our results are obtained via time-dependent nonequilibrium Green 
function. Both parallel and antiparallel lead magnetization alignments are considered. 
The main findings include a wave of spin accumulation and spin current that can 
change sign as the time evolves, spikes in the TMR signal and a TMR sign change due 
to an ultrafast switch from forward to reverse current in the emitter lead. 
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1. Introduction 

Spintronic[T] has proved to be of great technological importance with the development 
of memory storage devices and magnetic sensors based on giant magnetoresistance 
(GMR)[2] and tunnel magnetoresistance (TMR). [3] More recently, one of the major 
achievements in the spintronic field was the coherent control of single electron spins in 
quantum dots. 01 El El [TJ, [HI El EEH EE21 [13] Such a control is a fundamental step toward 
the generation of quantum bits based on the electron spin for further implementation 
of quantum computers. [HI [151 EH] There are a few recently developed techniques to 
manipulate coherently the electron spin in quantum dots. In those techniques square 
pulses of bias or gate voltages are applied in order to prepare and measure spin states. [17] 
This turns quite desirable the study of quantum transport in the presence of pulses of 
bias voltage. For single level quantum dots coupled to nonmagnetic leads, it is well 
known that coherent oscillations (ringing) of the current appear when a step like bias 
voltage is applied across the junction. [18] In the presence of ferromagnetic leads and/or 
Zeeman-split level, though, it was found that the ringing response of the current develops 
spin fingerprints. [T9l [20] Additionally, polarized current spikes can be generated when 
a bias voltage is abruptly turned off. [21] 

Systems with high TMR values are of current interest in spintronics, motivated 
mainly by magnetic memory and sensor applications. [22] Recently, it was found 
relatively high TMR values (80%-200%) in magnetic tunnel junction (MTJ) consisting 
of Fe and Co spaced by an oxide layer. [231 EU EH] In particular, if instead of an insulator 
layer we have a quantum dot sandwiched by two ferromagnetic leads, [26] it is possible 
to have additional effects, like spin-accumulation and strong Coulomb interaction that 
result in a more wealth physics of the TMR response. [2?J EH] [291 [301 ETJ [321 [Ml EH1 ^ 

In the present work we extend the recent study in Ref. [21] by considering a 
sequence of square pulses of bias voltage instead of only a single pulse. A study on 
electron spin dynamics when a sequence of voltage pulses is applied in a quantum dot 
system, was recently developed by Stefanucci. [37] Since square waves are of fundamental 
importance for the conventional electronics (e.g. digital switching circuits, synchronous 
logic circuits, binary logic devices) [38] we believe that it is relevant to consider the 
interplay between square waves of bias voltage and spin dependent phenomena. The 
calculation presented here is based on the nonequilibrium Green function technique. 
The main finding is a wave-like behavior of the spin accumulation, spin current and 
TMR, that can switch sign periodically as the time evolves. In particular, the TMR 
develops a periodic singularity that results in relatively high TMR and a fast switching 
of its sign. 

The paper is organized as follows. In Sec. II we derive an expression for the current, 
in Sec. Ill we show and discuss the results and in Sec. IV we conclude. 
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2. Hamiltonian and Transport Formulation 

To describe the system illustrated in Fig. [T]we use the following Hamiltonian 

H = e ka V (t)cl ari c kav + Y £d{t)dld a 

k<r»7 a 

+ J2( V( Lr, d " + V * d * C *°v), (!) 

where c^ ar) (c^ ) and d a (dj.) are the annihilation (creation) operators for electrons in 
the lead r\ and in the dot, respectively. The energies e^it) and edit) are the time- 
dependent energies for the electrons in lead rj (77 =L or R for left or right) and in the 
dot, respectively. The labels k and a denote the electron wave vector and the spin, 
respectively. More explicitly, these energies are written as e^ ar1 {t) = e£ + A v (t) and 
€d(t) = e d + A d (t), where e£ and e d are time-independent energies and A Vtd (t) gives 
the time evolution of the external bias. Using Eq. (CQ) inside the current definition 
I%(t) = —ie([H, NJ}]), where e is the electron charge (e > 0) and N% = ^2^^ c^ ari , we 
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Figure 1. Schematic energy diagram for the system considered. The dot level Cd{t) is 
coupled to a left (emitter) and to a right (collector) lead with chemical potentials 
fiL and fiR, respectively. During a pulse of bias voltage the level edit) and the 
right chemical potential \ir are shifted down, thus allowing charge and spin transport 
through the system. When the bias is turned off the level ed(t) becomes above both 
fiL and [ir, which results in a transient discharging of the dot. This discharge process 
is spin-dependent due to the ferromagnetic leads. 
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can show in the noninteracting[59"] case and wideband limit that [101 E] 

m) = -erHMt)) + J */,(e)M4n,(M)]}, (2) 

where f v (e) is the Fermi distribution function for lead r], and (n a (t)) is the time- 
dependent dot's occupation, given by 

(n a (t)) = lm{G< a (t,t)} 

= / |^)I^M)l 2 - (3) 

v J 

The function A ar] (e,t) is defined as 

A av (e,t)= f dhGUt^e^-^- 1 ^ 41 ^, (4) 

J — oo 

where the retarded Green function in the noninteracting model is given by 

G r aa (t, h) = -i6{t - tOe-^-^e"^! dIed{I \ (5) 

with T a = + and = e + Ad(t). The quantities and give the tunneling 
rate between left and the right leads into/out the dot, respectively. Substituting Eq. 
© into Eq. gD we find 



A av (e, t) = - i f dhe^+^-^e'^H di ^)-^ 



(*)] 



Solving Eq. (jBJ) for a bias voltage of the kind 

oo 

W*) = V v ° /d <j>(t) = V v ° /d J2^t - s n ^)6(s n -t), (7) 

n=l 

with A r? (t) = — V v and A^(t) = —Vd(t), we obtain 

4^(e,t) = -i f dt xl {tM)e- ia ^ m)+ ^=° { - ir+lsn] 



n=0 " s ™ 
f 



-i / ^(MiK^ - * 1 ^, (8) 

where iV gives the last instant tjy (time in which V(t) is turned on or off) before the 
time t in which A av (e,t) is being evaluated. The others quantities are defined as 

an = V$ ~ Vl (9) 
= e i{e+i ^- eo){t - tl \ (10) 
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and £ = for N = and £ = 1 for TV > 1. It is yet valid to mention that is assumed 
constant throughout the leads. Physically, this means that the electronic system has 
enough time to screen an external electric field as it evolves in time. This assumption 
implies bias modulations not faster than the typical plasma frequency. [42] The three 
integrals in Eq. (jSJ) can be solved analytically, thus resulting in the following expression 

e t(e- eo -H I f)t e -ia,W(t)+I^L (-l) B+1 «n] 
Arr;(e, t) = — p- h 

e - e + i^- 



N-l 



e i(e-e +i^)t e -i av it<j ) (t)+Y:fL n+1 (-iy +1 si] x 



n=0 

e - e + - c^/n 

1 _ p*[e-eo+*^-a>7^(t)](t-ajv) 

- - . (11) 

e - e + i^f - a v (fi(t) 

In this last equation f n = l for n even and for n odd. Using Eq. (fill) into Eqs. (J2]) and 
([3]) we determine the dynamics of the spin polarized current and the dot occupation, as 
described in Sec. 4. 



3. Parameters 



In our numerical calculations we have described the ferromagnetic leads via spin 
dependent tunneling rates T^, given by [34"] 

v L a = r [i + {-1)^VlI (12) 

= r [i ± (-i)^M (is) 

where To is the lead-dot coupling strength, pi and pr are the left and right lead 
polarization, and the + (— ) sign in corresponds to parallel (antiparallel) magnetic 
alignment of the leads. In particular, in the present work we assume Pl = Pr = 0.4. 
The others quantities involved are the dot's level eo = 5Po, the voltages intensities 
Vl = 0, Vr = 20r and = 10r , the pulse width and the inverval between pulses, 
s n — s n _i = 3h/T , and the temperature ksT = T . In what follows we present the 
results. 



4. Results 



4-1. Spin accumulation 

Figures [2]^a)-(b) show the spin resolved dot occupations in both (a) P and (b) AP 
alignments. When the bias is turned on (t = 0, 6h/T , 12/i/T , ...) the occupations 
increase in time due to the resonant condition fiR < £d{t) < ^(t) that is achieved 
inside the pulse length. In the parallel configuration presents a stepper enhancement 
compared to n±. This is related to the inequality > T^, that gives a faster response 
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Figure 2. Dynamical spin-resolved electronic populations in the (a) parallel and (b) 
antiparallcl configurations and (c) the spin accumulation m = n-\ — n^. While in the 
parallel configuration both n-f and attain equal plateaus in the stationary limit, in 
the antiparallcl case we find > in this limit, which results in a higher amplitude 
of the spin accumulation wave in the AP than in the P alignment. Additionally, the 
spin accumulation wave changes sign in the P case. 

for the incoming up spins. As the time evolves n-f and nj, tend to the same value, inside 
the pulse. In contrast, in the AP alignment nj tends to saturate above n±. These 
asymptotic behaviors can be easily understood in terms of the tunneling rates in both 
alignments. While in the P case = T^, thus resulting in n-f = in the stationary 
limit, in the AP configuration > T? and rf < Ff, which gives rise to > in 
the stationary limit. When the bias voltage is turned off (t = 3h/T , 9h/T , 15h/T , ...) 
the level becomes above /i£ and /iR, thus the spin populations in the dot discharge 
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Figure 3. Spin current agains time in both parallel and antiparallcl alignments and in 
both ferromagnetic leads. In the left lead (emitter) lf pin develops higher values than 
I^ p F in inside the bias pulse. In contrast, when the bias voltage is turned off lf p f n attains 
much higher negative values than lf pin . 



toward the leads. The minimum values achieved by nj and during the discharge 
process come from thermal excitation. Note that just before a pulse of bias voltage we 
find n-f > nj in the P case while n-f = in the AP configuration. This is related to the 
broadening of the dot level that is spin dependent, even though it is spin degenerate. 
Since in the P case + > rf + the level broadening for spin up is larger than 
for spin down. This makes thermal excitations more pronounced for spins up than for 
spins down, thus resulting in an equilibrium spin accumulation. In the AP alignment 
since r * + = Tf + the thermal excitation is equally distributed for both spins. In 
Fig. [2](c) we show the spin accumulation in the dot, m = rij —ny. For AP configuration 
m changes intensity preserving its sign. In contrast, in the P alignment m oscillates 
between positive and negative values as the time evolves. The sign reversion of m in the 
P case comes from the faster charge/discharge of spin up electrons in the dot, compared 
to the spin down electrons. 

4-2. Spin current 

Figure^a) shows the spin currents (I sp in — h~ -Tl)@3] m the emitter (left) and collector 
(right) lead in both magnetic alignments. Observe that just after a bias voltage is turned 
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Figure 4. Tunnel magnetoresistance against time in the left (solid line) and right 
(dotted line) leads. In the time range in which V(t) = the TMR increases up to 
200-300% due to the relatively slow discharging of spin down electrons in the dot in 
the parallel alignment. This sustains I p much longer than I , thus enhancing the 
TMR in time. Additionally, when the bias voltage is turned off, a singularity in the 
TMR is observed. This is related to the zero value attained by the current in the left 
lead when it is passing from forward (positive) to reverse (negative) current. 



on, both I p pin and I Ap n are on top of each other. As the time evolves I Ap n is suppressed 
to zero while I p vin tends to a nonzero stationary value. In contrast, when the bias 



voltage is turned off, I AP n assumes transient values higher in modulos than I p pin . The 
negative values of the spin currents just after a pulse end means that the spin current 
is discharging into the emitter lead with a majority up component. [44] In contrast, in 
the collector lead [Fig. [31(b)] I^f n is aproximatelly zero throughout time, with small 
oscillations (ringing like) whenever a bias voltage is turned on or off. In contrast, I p pin 
assumes relatively high negative values. Comparing Figs. EJ^a) and [3](b) we note that 
in a transient regime the spin current that leaves the emitter does not arrive in the 
collector lead. Similarly to the total current (J = /j + that has a continuity equation 
given by 

lL + lR = 6 % (14) 
where n — + the total spin current satisfies the following continuity equation 

tL tR dm , . 

^ spin ' ^ spin e ^ i (l^/ 

which was used to check the accuracy of our numerical results. 



4.3. TMR 

Figure H] shows the TMR= (I p — I AP )/I AP against time in the left (solid line) and right 
(dotted line) leads. Here I p and I AP means the total current (7| + 1±) in the parallel 
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Figure 5. TMR and total currents I p and P 4 ^ 3 against time around t = 3h/To- The 
TMR changes sign when I AP switches from forward (positive) to reverse (negative) 
current. 



and antiparallel configurations, respectively. In the time range in which V(t) — the 
TMR increases up to 200-300% due to the relatively slow discharging of spin down 
electrons from the dot into the leads in the parallel configuration. This sustains I p 
much longer than I AP , thus enhancing the TMR along the time. Note that in the 
parallel case the spin down electron in the dot is weakly coupled to both leads (majority 
spin up population in both sides), while in the antiparallel case both spin components 
are strongly coupled to at least one electrode. This turns the spin down discharging 
process slower in the parallel case, thus resulting in a maintenance of the total current 
I p for longer times compared to I AP '. It is valid to mention that whenever a bias voltage 
is turned off the TMR develops a singularity as described in the next figure. 

Figure [5] shows the TMR and the total currents I p and I AP in a time range around 
3H/Tq. The divergence and sign change of the TMR observed in this figure is related 
to the spin-dependent discharging process of the dot. Note that after the bias voltage 
is turned off (at 3h/T ) the currents I p and I AP transiently pass from direct (positive) 
to reverse (negative) currents. During this transition the current I AP crosses the zero 
value before I p and attains a slightly more negative value than I p . This turns into a 
singularity in the TMR and a change of its sign. 

5. Conclusion 

We calculate spin polarized current, spin accumulation and TMR in a quantum dot 
coupled to ferromagnetic leads in the presence of a bias voltage that evolves in time 
as a square-wave. We report a wave like spin accumulation that switches sign in each 
period when the leads are parallel aligned. We also found a much larger spin current 
in the antiparallel configuration during the transient discharging process, which quite 
contrasts with the stationary regime where we find an antiparallel spin current equal 
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zero. Finally, we find a large enhancement of the TMR in the time range in between 
bias voltage pulses. This is due to the slow discharge of the spin down electrons in the 
parallel configuration. A sign reversion of the TMR is also observed in the emitter lead, 
which is related to a sign change of the current. 

The authors acknowledge I. Larkin, R. Zelenovsky, and A. P. Jauho for valuable 
comments. This work was supported by the Brazilian Ministry of Science and 
Technology and IBEM (Brazil). 
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